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of  the  minimum  model  size  necessary  for  the  existence  of  a  turbulent 
wake  at  collapse. 
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Abstract 


The  modeling  of  the  turbulent  wake  of  a  self-propelled  body  in  a 
stratified  fluid  is  discussed.  The  scaling  parameter  is  taken  to  be  the  internal 
Froude  number,  using  the  speed  and  diameter  of  the  body  and  Vaisala  frequency 
of  the  fluid.  Arguments  are  presented  to  suggest  that  the  ratio  of  the  time  for 
the  wake  to  collapse  to  the  characteristic  time  of  the  turbulence  is  invariant 
with  model  scale  if  the  wake  is  turbulent  up  to  and  including  the  collapse.  A 
number  of  criteria  for  the  existence  of  a  turbulent  wake  are  discussed.  Numeri¬ 
cal  estimates  are  made,  assuming  typical  values  of  the  model  speed  and  of  the 
Vaisala  frequency,  of  the  minimum  model  size  necessary  for  the  existence  of 
a  turbulent  wake  at  collapse. 
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Model  Scaling 


In  this  note,  the  modeling  of  the  turbulent  wake  of  a  self-propelled 
body  of  length  L_  traveling  at  a  speed  U  in  a  stratified  fluid  which  is 
characterized  by  a  Vaisala  frequency  f\f  is  discussed.  The  diameter  of  the 
body  is  0  ,  the  diameter  of  the  propeller  is  c/  and  it  will  be  assumed 

that  L/D  =  JO  and  [)/ c/=  2  so  that  L/J*ZO.  For  the  full  scale  or  proto¬ 
type  body,  typical  numerical  values  are  =  6  knots  3m/ sec,  L.-  /CO 
and  A/  =  lO’^sec-^. 

The  scaling  parameter  is  the  Froude  number 

F  =  U/fs/O) 

which  is,  of  course,  equivalent  to  Richardson  number  scaling  since 

Ri  =  >/Fz . 

In  the  usual  way,  U  .  0  and  A/  are  chosen  so  that 


where  the  subscripts  p  and  A'/  refer  to  the  values  of  the  parameters  for 
the  prototype  and  model  respectively. 

In  both  the  model  and  prototype,  the  physics  must  be  the  same  for  a 
meaningful  model  experiment.  The  propeller  produces  turbulence.  The  wake 
grows  by  entraining  fluid  with  a  consequent  decrease  in  the  turbulent  energy 
and  increase  in  the  potential  energy  of  the  fluid  within  the  wake.  Simultaneously 
there  is  viscous  dissipation  of  the  turbulence  energy.  Finally,  when  the  turbu¬ 
lent  energy  drops  to  a  sufficiently  low  level,  the  wake  collapses. 
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One  of  the  consequences  of  Ko's  model  of  the  growing  and  collapsing 

* 

wake  (Ko,  1971)  is  that  the  time  for  the  wake  to  collapse,  tc  ,  is  nearly 
independent  of  the  Froude  number  and  is 

tz  =  z/a/. 

This  is,  of  course,  in  agreement  with  the  available  experimental  data 
(Schooley  and  Stewart,  1963;  Van der  Watering,  1966).  Another  important 
time  scale  is  that  of  the  turbulence  (Batchelor,  1967) 

tT~$/U 

where  is  the  integral  scale  of  the  turbulence  and  CA.  is  the  turbulent 
intensity. 

Clearly,  the  ratio  between  these  times  must  be  equal  for  both  model 
and  prototype  if  the  phenomenon  of  wake  growth  and  collapse  is  to  be  cor¬ 
rectly  modeled.  This  ratio  is 

tz/tT  =  2U./(MJ) 

The  turbulent  intensity  is  a  fraction  of  [J  , 

U  =  *  U  J  G  <  °<  *  /  . 

Ko  suggests  that  0(  ^  0<  2  £"  just  behind  the  propeller  and  this  is  in  rough 
agreement  with  the  results  obtained  by  Naudascher  (1965).  It  is  also  to  be 
expected  that  the  integral  scale  is  a  fraction  of  the  propeller  diameter,  qJ  , 

i.  e.  , 

J-f  J  ,  OCfCl. 


-3- 


Therefore 


That  is,  the  ratio  of  the  time  for  the  wake  to  collapse  to  the  character¬ 
istic  time  of  the  turbulence  will  have  the  same  value  for  both  the  model  and  the 
prototype  if  the  turbulent  intensity  is  the  same  fraction  of  the  body  speed  and 
the  integral  scale  is  the  same  fraction  of  the  propeller  diameter  in  both  model 
and  prototype. 

In  both  model  and  prototype  the  self -propulsion  is  accomplished  by  a 
propeller;  both  the  model  and  prototype  propellers  are  geometrically  similar 
and  the  ratio  of  the  rotational  speed  of  the  propeller  to  body  speed  is  the  same 
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for  both  model  and  prototype.  Since  both  the  geometry  and  the  velocity 
triangles  of  the  blades  are  similar,  it  is  expected,  provided  that  the  model 
propeller  produces  a  turbulent  flow,  that  the  ratios  of  the  integral  scale  to 
the  diameter  and  the  turbulent  intensity  to  the  body  speed  are  invariant  with 
scale.  That  is,  it  can  be  expected  that  0(  and  &  are  truly  constants. 


Then  the  ratio 
mated  if  it  is  assumed  that 


also  invariant.  Finally, 
0(  =  1/4  and  (3  1/2, 


this  ratio  can  be  esti- 


/  oo  . 


These  estimates  suggest  that  the  collapsing  wake  can  be  properly 
modeled  at  any  scale  provided  that  the  wake  is  turbulent  at  that  scale.  But, 
is  the  wake  turbulent  behind  the  propeller  and,  say,  up  to  the  point  of  collapse? 
It  is  clear  that  the  Reynolds  number  will  be  very  different  in  model  and  proto¬ 
type,  but  there  must  be  a  similarity  in  that  a  buoyancy  range  and  an  inertial 
range  must  be  present  in  the  turbulent  spectrum  at  model  scale  if  they  are 
present  at  prototype  scale. 

Therefore  consider  the  conditions  at  two  sections;  just  behind  the 
propeller  and  at  the  point  of  collapse.  Taking  the  time  of  collapse  as 

tc  =  2-/V, 

the  point  of  collapse  is 

Xc  =  -  2.  U/A/. 


-5- 


Therefore 


using  the  values  of  U  ,  D  and  /V  for  the  prototype.  So  the  state  of 


the  flow  must  be  estimated  at 


K  -  o 


it  =  /ZO  d 


In  all  the  numerical  work  below,  it  will  be  assumed  that  the  speed  of  the 


model 


,  U„  =  1  m/sec. 


Conditions  for  the  Existence  of  a  Turbulent  Wake 


A  necessary  condition  for  the  existence  of  an  inertial  range  in  the 
turbulent  spectrum  (Phillips,  1966)  is  that 

n'h  f  u*  )'/*  , 

R,  =(ir)  >>f- 


It  will  be  assumed  that,  just  behind  the  propeller  (  X  =  0), 

u  =  U/¥  j  J  =  «//z, 

then  /. 

/£-(¥)*(&)  YfA 

With  (J  =  1  m/sec,  J/L  -  1/20,  and  2^  =  10“^cm^/  sec, 

Rf  *  9  L'x 
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with  L  in  cm.  It  is  clear  that  any  model  longer  than  about  50  cm  easily 

P  'A 

satisfies  She  condition  K/0  »/. 

Next,  estimate  f^J  at  X  =  120  c/  .  Estimates  of  U  and 
at  7t  =  120  d  arc  required.  The  only  experimental  results  which  are 
available,  those  of  Naudaschcr  (1965),  will  be  used.  Using  the  data  in 
Naudascher's  Table  3  and  Figure  15 


Jf/J  *  o.s  j  u/  u  z  7  k  /o 


•3 


Therefore 


/?. 


-  ( 3 .  3  </»  ~J  (  y/  (  L  / 


/Zo 


Z  /.3  L 


>/x 


with  in  cm.  Thus,  for  the  existence  of  an  inertial  range, 


L  ‘h  »  /, 


and  it  appears  that  L*  -  100  cm  will  be  adequate.  Therefore  there  will  be 
an  inertial  range  at  collapse  if  L.  ^  1  meter. 

Next,  under  what  conditions  will  a  buoyancy  range  exist?  It  is  clear 
that,  just  as  in  considering  ft*  ,  the  crucial  conditions  occur  at  %  =  120  d 
The  buoyancy  wavenumber  Kt  which  separates  the  buoyancy  range  and  the 


*  Naudaschcr  uses  O  for  the  diameter  of  the  turbulence  generator 
and  L.  for  the  integral  scale;  our  d  and  £  ,  respectively. 
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equilibrium  range  is  (Lumlcy,  1964;  Phillips,  1966) 

K  «  ("’A)"'' 

where  f  is  the  rate  of  energy  dissipation.  For  a  buoyancy  range  to  exist 

( 'A?  )  » 

where  . 

/  >  /  \ 

J  =  (  */*) 

is  the  Kolmogorov  micro-scale.  Taking 

e  = 

this  condition  is 

u3/(yJ/ Vz)  »  /. 

This  can  be  written  as 

mmxtxv)  rx » '. 

or  with  LA.  f  U  =  7  x  10  ^ ,  I/d  =  0.  5  from  Naudascher's  results,  and 
-  1  m/  sec. 

L  >>/ 

with  Z«  in  cm.  Again,  it  seems  reasonable  that  Z.  ^  1  meter  will  be 
*)< 

adequate.  *“ 


*  At  X  “120  d  ,  the  ratio  is  actually  larger  than  ^  by  the  amount 

of  energy  extracted  in  buoyancy  range.  Including  this  consideration,  this 
condition  is  given  bv 

L  »>  /. 

Again,  it  seems  reasonable  that  L,  100  cm  is  adequate. 
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Finally,  it  might  be  required  that  the  Reynolds  number  based  on  the 
dissipation  length  scale,  ^  , 

Again  consider  only  the  case  X  =  120  J  .  Naudascher  (Figure  13)  gives 
curves  of  X/xw(  X  c«o  /*-<>•  135)  as  a  function  of  V  ( -  X0J ^ U 
Assuming  that  *0=o, 


with  in  cm.  For  / £  /.  £■  / O ^  ^Jies  in  the  range  0.  5 
to  0,  2.  Therefore  it  will  be  assumed  that 

X/X«o  -  C>.¥  and  \~/</  *  O.ZJS- 


to  obtain  an  estimate  of  .  Therefore, 

A 


a,  -  mt) 

=  o./&>  L 


i-vjyuL 


d  L  n  y 


with  ^  in  cm.  Again,  it  appears  that  L*  1  m  will  be  satisfactory. 
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Conclusions 


a)  The  ratio  of  the  time  for  the  wake  to  collapse  to  the 
time  scale  of  the  turbulence  (  /tT)  wi  n  scale 
correctly  with  the  Froude  number  F"-U/t^O  if  the 
wake  flow  is  turbulent  at  model  scale. 

b)  For  the  model  wake  to  be  similar  to  that  of  the 
prototype,  the  requirements  are: 

for  an  inertial  range, 

R'f*  [(“/v)  (*/J)  ( JA)  /?4  ]  »  /; 

for  a  buoyancy  range, 

(^/u)3(  d/s)  (  o/jXd/l  )  r  zRl  » tj 

and  finally, 

/?>  =  ( “/u)( d/L )Rl  »/ 

with 

Rl  -  UL/v . 

Using  Naudascher's  results,  in  all  cases  these  inequalities  can  be 
satisfied  by  taking  RL  =  106,  i.  e.  ,  (J  =  lm/ sec.  L—  =  1  m  with 
"J/  =  10  ^cm^/sec. 
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